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o : 

Abstract. This work applies a Rayleigh-Brillouin light spectrum analysis in order to establish 
q 1 a causality test by means of a frequency cone. This technique allows to identify forbidden and 

qj | unforbidden regions in light scattering experiments and establishes if a set of linearized transport 

. equations admits causal solutions. It is shown that, when studying a relativistic fluid with its acoustic 

modes interacting with light, Eckart's formalism yields a non causal behavior. In this case the 
solutions describing temperature, density and pressure fluctuations are located outside the frequency 
cone. In contrast, the set of equations that arises from modified Eckart's theory (based on relativistic 
q ' kinetic theory) yields solutions that lie within the cone, so that they are causal. 
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^ : 1. INTRODUCTION 

*n : 

It is well-known that Eckart's formalism contains generic instabilities while describing 
linear fluctuations in a relativistic fluid [1]. This fact is due to the acceleration term 
present in the constitutive equation for the heat flux J' q [2]. 
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In this paper we focus on the causality issue with a novel approach. Since light inter- 
acts with traveling waves in a fluid, absorbtion and emission processes cause frequency 
shifts depending on the wave speeds. The possible values of these shifts are restricted by 
special relativity and its description depends on linearized relativistic hydrodynamics. 
In this paper we compare the predictions of the intensity spectrum of scattered light off 
acoustic modes as given by Eckart's formalism and the modified relativistic hydrody- 
namics as included in Ref. [3]. 

In section II the balance equations for a relativistic fluid are linearized using Eckart's 
formalism, and the corresponding spectrum equation is calculated. In section III ana- 
log calculations are performed using the modified Eckart's formalism. In section IV the 
Rayleigh-Brillouin spectrum is described as the test methodology is proposed, introduc- 
ing the light-frequency cone for fluid perturbations. This section also includes compar- 
isons between the spectra arising from sections II and III. Section V is devoted to final 
remarks. 
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2. LINEARIZED TRANSPORT EQUATIONS WITH ECKART'S 

FORMALISM 



For a simple, relativistic, non-degenerate fluid we may work with the balance equations 
as given below [4] [5] [6]. 
The continuity equation: 

-|» + n0 = O, (2) 
at 

where is defined, as usual, as uY v = 0, and n is the number of particles per unit of 
volume. 

The momentum balance equation: 

(9 jc" i9T kT <9 2 A 

p-(0) +nkV 2 (T)+kTV 2 (n) -AV 2 (0) - ? V 2 — - = 0, (3) 

where /: is the Boltzmann constant, K is the thermal conductivity and A is the longitudinal 
kinematic viscosity [7]. The quantity p = jj(ne + p) is identified with the local enthalpy 
of the fluid. Use has been made of the equation of state p = nkT. 
And, finally, the energy conservation equation: 

with Dj the thermal diffusivity [7]. 

For a proper linearization of this set of equations we decompose each thermodynamic 
function (T, n and 9) in terms of its value at equilibrium and a fluctuation term [5], that 
is 

X = X + 8X (5) 

Subscripts "o" refer to values at equilibrium. Making this substitution and neglect- 
ing second order terms in fluctuations, the following set of linearized equations, with 
Eckart's formalism arises. For particle conservation: 

j- t (8n)+n o 80 =0 (6) 

For the momentum balance: 

fiJW +k TVHSn) + nkV\5T) - A*(M) - - = « P) 

And the conservation of energy equation, making the substitution D T = m * c , reads: 

| (5T) _?Zk 59 _^v^)-^^ = 0, (8) 
at 3 mn () c n mn„c„c z at 

where c n is the specific heat at constant particle density. 



Light Cone analysis of relativistic first-order in the gradients hydrodynamics December 3, 2010 2 



In order to obtain the spectrum [7], Fourier and Laplace transforms are performed 
with the set of equations (6-8). Then, the following dispersion equation arises (for further 
details see reference [3]): 
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The evaluation of the determinant leads to a quartic polynomial that reads: 

d 4 s 4 + s 3 + d 2 s 2 q 2 + s(d 3 q 4 + d Q q 2 ) + d^q 4 = 0, 
the coefficients d\, d-$, di, d\ and do are given by 
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where c p = yc n is the specific heat at constant pressure. 

The reader may recall that the fourth root of Eq. (10) doesn't allow a finite spec- 
trum, since the corresponding Fourier transform is divergent [8]. Physically, Eckart's 
framework predicts a significant (unobserved) broadening of the central peak, directly 
associated to the thermal conductivity coefficient present in Eq. (11). This coefficient 
appears due to the heat acceleration coupling, Eq. (1) . 



3. FLUCTUATION TRANSPORT EQUATIONS WITH MODIFIED 

ECKART'S FORMALISM 

If Eq. (1) is replaced in the basic formalism with the alternative constitutive equation, 
directly obtained from kinetic theory [3]: 

Jq = —LttT* —Lnjn* (16) 

The linearized continuity equation (6) remains unaltered, while the linearized momen- 
tum balance equation is now [3]: 

p /jm +kTvHSn) + nkvHST) _ A ^ M) _!*T^jm_ % w = o 

(17) 
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and the linearized energy balance reads: 

4-(5r) - ^80 - ^-V 2 (8T) - ^-V 2 (8n) = 
at 3 n c n n a c n 

Ljt is an effective thermal conductivity given by 

PTo , 



L TT = k(1 + 



c 2 p k T 



where kj is the thermal compressibility and j5 the thermal-expansion coefficient. 
Defining z = in the case of an ideal gas, Eq. (19) reduces to 

L tt = k(1+z) 



(18) 



(19) 



(20) 



for z « 1. The other coefficient, L n r, has no counterpart in the classical regime and is 
given by 



kT„ 



n k T c 2 p 



(21) 



For z<<1, equation (21) reduces to 



L nT =Z-K 



(22) 



for an ideal gas. As in the previous section, Fourier and Laplace transforms are taken in 
the whole system, so that the following dispersion equation yields: 



s no 
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In this case we don't have a quartic, but a cubic polynomial equation: 

s 3 + b 2 s 2 q 2 + s(b 3 q 4 + biq 2 ) + b 4 q 4 = 0, 
with the coefficients given by: 

A L T t m T 2 n () n L nT 



p n c„ c^po 
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(23) 



(24) 

(25) 

(26) 
(27) 
(28) 
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The Rayleigh peak's width [3] reads now: 

Ar = T^ = — ( Ltt 1 — )' ( 29) 

b\ cpn () k 

In this last expression, the pathological broadening that arises with Eckart's formalism 
is absent. 



4. THE RAYLEIGH-BRILLOUIN SPECTRUM AND THE LIGHT 
FREQUENCY CONE FOR FLUID PERTURBATIONS 

The reader may recall that in a simple relativistic fluid two transverse and three lon- 
gitudinal hydrodynamics modes are found [7] [9]. The latter are at the origin of light 
scattering, giving rise to what it is called the Rayleigh-Brillouin spectrum, which con- 
sists of a central peak (the Rayleigh peak) and two secondary peaks (the Brillouin peaks) 
that appear from a Doppler effect. Denoting (Oq the frequency of the center of Rayleigh 
peak, then the two Brillouin maxima are located at COq ± c s q, where c s is the speed of 
sound and q is the wave number. 

When waves are present in a fluid, the causality properties can be verified by making 
it interact with light. If light with a frequency v is re-emitted with a frequency v' = 
v(l + (where u is the fluctuation's propagation speed), then v' would never become 
greater than 2v. Thus, a light cone on the frequency domain is generated. 

The di\ coefficient is the clue to determine that Eckart's formalism yields a non-causal 
solution [8] [10], since the width of Rayleigh' s peak is now enhanced by the amount: 

A ^ (30) 
KT 

The term cq delimits the frequency cone so if the following relation holds 

P° cA ^ ri-\\ 
^f >Cq (31) 

then the solution is out of the cone and, thus, is non-causal. 

Following the same reasoning included in Hiscock and Lindblom's work [1] some 
calculations for water at a temperature of 293 K and pressure of 1 bar are shown below: 

^ ~ 10 34 — (32) 
kT sec 

evaluating cq with the typical value q = XQpm^ 1 

1 

sec 

Since 10 34 >> 10 15 the solution is clearly out of the frequency cone, and thus is non- 
causal. Physically, this would mean that the interacting particles of the fluid move with 
a supraluminical speed after collisional processes. 



C4-10 15 — (33) 
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5. FINAL REMARKS 



A test based on light scattering experiments has been developed to determine if a set 
of linearized hydrodynamic transport equations yields to a causal solution or not. When 
performing the calculations under Eckart's formalism, and increasing z, the Rayleigh 
peak's width invade the two Brillouin ones, entering the forbidden region of the light 
cone. It is clear that the non-causality issue arises when coupling the heat flux with the 
acceleration. This work has been developed with first order in the gradients theories. 
With the standard version of extended theories [11], the fourth root (s 4 ) is conserved, 
but damping terms are included. It is not clear if this damping preserves the overlapping 
of the Brillouin peaks mentioned above. Also, the damping coefficients cannot be kinet- 
ically supported by a Chapman-Enskog expansion. These issues deserve a closer look 
and will be studied in the near future. 

The authors wish to thank A.L. Garcia-Perciante and L.S. Garcia-Colin for their 
valuable and encouraging comments. 
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